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Abstract 

A connected algebraic group in characteristic is uniquely determined by its Lie 
algebra. In this paper an algorithm is given for constructing an algebraic group in 
characteristic 0, given its Lie algebra. Using this an algorithm is presented for finding 
a maximal reductive subgroup and the unipotent radical of an algebraic group. 

1 Introduction 

Due to a wide variety of applications of algebraic groups there is considerable interest in 
doing computations with such groups. For example, [3], [Tj contain both applications of 
algebraic groups and computational methods. The most natural way to give an algebraic 
group is by a set of polynomial equations. However, it is not easy to determine properties of 
the group from these equations (for example, whether it is simple, solvable, or nilpotent). 
A standard way to get around this is to study the Lie algebra of the group. This Lie 
algebra reflects many of the properties of the group, and can therefore be used to get a 
picture of the structure of the group. However, by looking at the Lie algebra alone it is for 
instance not possible to construct subgroups, such as the unipotent radical, or a maximal 
reductive subgroup. It is the aim of this paper to provide methods for constructing such 
subgroups by taking the inverse route. In other words, we describe algorithms for going 
back to the group from its Lie algebra. 

For this a few restrictions are necessary. Firstly, in characteristic different from zero 
it is in general not true that to a given Lie algebra corresponds a unique algebraic group. 
There can be more than one. So we restrict to fields of characteristic zero. Second, since 
the Lie algebra of an algebraic group is equal to the Lie algebra of the connected component 
of the identity, we restrict to connected algebraic groups. We note that, also with these 
restrictions, not every Lie algebra corresponds to an algebraic group. Lie algebras that do 
are called algebraic. There are methods to decide whether or not a given Lie algebra is 
algebraic (cf. [1]). So we assume that we have given an algebraic Lie algebra, and we want 
to construct the corresponding group. 
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Now we state the problem more precisely. Let F be a field and let F be its algebraic 
closure. By R(n, F) we denote the polynomial ring F[xn, . . . , x nn ] (n 2 in determinates) . A 
subgroup G C GL(n, F) is said to be algebraic if there is a set of polynomials S C R(n, F) 
such that G consists of the g G GL(n, F) with f(g) = for all f E S. Let G be an algebraic 
group, and let 1(G) C i?(n, F) be the ideal of all polynomials vanishing on G. Then G is 
said to be defined over F if 1(G) is generated by polynomials in R(n, F). 

In this paper we take F to be a field of characteristic (Q for example), and consider 
algebraic groups that are defined over F. Since F is perfect, any subgroup C GL(n, F) 
defined as the vanishing set of a set of polynomials in R(n, F) is automatically defined over 
F (cf. [8], §34 ). We call a set S C R(n, F) such that G consists of all g where all elements 
of S vanish, a set of defining polynomials for G. (We do not assume that they generate 
1(G).) 

The Lie algebra g = Lie(G) of an algebraic group G C GL(n, F) is a subalgebra of 
gi(n, F) (the Lie algebra of all n x n-matrices with coefficients in F). Now if G is defined 
over F, then g is a subalgebra of gl(n, F). So the problem that we consider is the following. 
Given an algebraic Lie subalgebra g C gl(n, F), compute a set of defining polynomials (over 
F) for the unique connected subgroup G C GL(n, F) such that g = Lie(G). 

The algorithm outlined in this paper consists of a few subalgorithms. In Section [2] we 
describe an algorithm that constructs the smallest algebraic group containing two given 
algebraic subgroups Gi, G 2 C GL(n, F). In Section [3] an algorithm is given for constructing 
the algebraic group corresponding to a Lie algebra that consists of nilpotent matrices. Let 
X G g; then by G(X) we denote the smallest algebraic subgroup of G such that its Lie 
algebra contains X. In Section H] we give an algorithm for constructing G(X) in case X is 
semisimple. Now these algorithms together solve the problem. Indeed, since g is algebraic, 
it is closed under Jordan decomposition. So g has a basis consisting of elements that 
are either nilpotent or semisimple. So for every basis element X we can construct G(X). 
Let X, Y be two basis elements. Then with the algorithm of Section [2] we can construct 
the smallest algebraic group H that contains G(X) and G(Y). The Lie algebra of H is 
generated by X, Y. Then we take a third basis element Z, outside of this subalgebra, 
and form the smallest algebraic group that contains H and G(Z). Continuing like this 
eventually we find G. 

In Section [5] we give an algorithm to decompose the Lie algebra g into a direct sum 
of a reductive subalgebra, and an ideal consisting of nilpotent elements. Together with 
the algorithms of the previous sections this yields an algorithm for finding the analogous 
decomposition of an algebraic group. This is one of the steps needed in the algorithms of 
[7]. Finally, the last section describes experiences with an implementation of the algorithm 
in the computer algebra system Magma. 

2 The algebraic group generated by two subgroups 

Let G%,G 2 C GL(n, F) be two algebraic groups defined over F. Let G C GL(n, F) be the 
smallest algebraic group containing both G\ and G 2 - Then also G is defined over F (this 
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follows for example from the algorithm given below). We consider the problem of finding 
a set of defining polynomials in R(n, F) for G. 

First we note that the multiplication map ■ : GL(n, F) x GL(n, F) — > GL(n, F) is 
a morphism of algebraic varieties. Using Grobner basis techniques we can compute the 
Zariski closure of the image of this morphism (see, e.g., [6], §1.8.3). This algorithm only 
uses operations in the base field. Hence the closure of the image of a variety defined over 
F is defined over F as well. So we can compute defining polynomials in R(n, F) for the 
Zariski closure G1G2 of the set G1G2 = {giQ2 \ 9% G G%\- Now from j2] (Chapter II, §7, 
Corollary 3), we have the following theorem. 

Theorem 1 (Chevalley) Let G\, . . . ,G S be connected algebraic subgroups 0/ GL(n, F). 
Then the group H generated by G\,...,G S is algebraic and connected. Moreover, there 
exists an m > such that every element of H can be written as a product of m elements, 
each belonging to aG{. 

Proposition 2 Let G±, G2 be connected algebraic subgroups o/GL(n, F). Let H C GL(n, F) 
be the group generated by G\,G2- Let A C H and suppose that 1 G A and that AG\ = A 
and AG~2~ = A. Then A = H. 

Proof. Let a G A, and g G G\. Then ag G AGi, and consequently ag G A. In the same 
way we see that ag G A if g G Gi- Let h G H. Then by Theorem [U h can be written as 
a product of elements from Gi U G 2 . It follows that ah & A for all a G A and h G H . So 
since 1 G A we conclude that A = H. □ 
This leads to the following algorithm for computing G: 

Algorithm 3 Input: algebraic subgroups Gi,G2 C GL(n,F). 

Output: defining polynomials for the smallest algebraic subgroup of GL(n, F) containing 

Gi, (?2- 

1. Let G be the trivial subgroup of GL(n, F). 



2. 


SetG' 


:= G. 


3. 


SetG' 


:= G'G X 


4- 


Set G' 


:= G'G 2 



5. If G' 7^ G then set G := G' and return to 2. Otherwise return G. 

Proof. We show that the algorithm is correct and that it terminates. The identity is 
always contained in G' . Secondly, if in a round of the iteration G' does not change then 
G'G\ — G' = G'G2- So by Proposition [2j G' = H, where H is the group generated by 
G\, G%. 

After k rounds of the iteration G1G2 • ■ • G1G2 (k factors G1G2) is a subset of G. So by 
Theorem [1] there is an m > such that after m rounds of the iteration we have G = H. But 
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then G' will not change in the next round of the iteration, and the algorithm terminates. 
□ 

3 The nilpotent case 

Let N C Ql(n,F) be a matrix Lie algebra consisting of nilpotent matrices. Then N is an 
algebraic Lie algebra ([]], Chapter V, Proposition 14). By G(N) C GL(n, F) we denote 
the connected algebraic group having Lie algebra N. We consider the problem of finding 
defining polynomials for G(N). 

Let N\, . . . , N r be a basis of N. Consider the polynomial ring P = F[T\, . . . , T r , Xij \ 
1 < i, 3 < n>]- Form the matrix M = exp(TiiVi + • • ■ + T r N r ), and let / be the ideal of P 
generated by X{j — M(i, j). By elimination techniques using Grobner bases we can compute 
the ideal J — I D F[xij}. 

Proposition 4 Set J' = {/ G Ffaj] \ f(a) = for all a G G(N)}. Then J = J'. 

Proof. First of all we note that every element of G(N) can uniquely be written as 
exp(J2i a iNi), where on G F ([lj, Chapter V, Proposition 14). Now let / G J and 
a = exp(^2 i aiN i ) G G(N). We show that f(a) = 0. Consider the ring homomorphism 
(p : F[Ti, . . . ,T r , Xy] — > F[xij], defined by sending Tj to «j, and leaving the x^ alone. Then 
(p(h) = h for all h G J. Also (p(h)(a) = for all h G /, as this holds for the generators of 
I. Hence f(a) = <p(f)(a) = 0. It follows that J C J'. 

Now let / G J', then /(exp(TiJVi + • • • + T r iV r )) = because f(J2i a i N i) = where 
the oti are arbitrary elements of F. Since x^- = M(i,j) + (x^ — M(i,j)), every polynomial 
in -Ffxjj] can be written as the sum of a polynomial in Ti, . . . , T r and an element of /. In 
particular, / = p(Ti, . . . , T r ) + h for a h G J. Then we substitute x^- i— >■ M(i,j). It follows 
that p = and / G J. Hence / G J. □ 
Remark. One of the advantages of this approach is that it returns the vanishing ideal 
of G(N). In [7J a more direct method is given. However, it does not have the same 
advantage. It works as follows. Let pi G R(n, F) for 1 < i < in be the linear polynomials 
that define iV as a linear subspace of gl(n,F). Let X = (x^j) be the matrix consisting of 
the indeterminates x^. Then G(N) is defined by the polynomial equations 

(X - l) n = and ^(log*(X)) for 1 < i < m, 

where log*(X) = — X^] 1 ~ X)*. However, the resulting polynomials are all of degree 
n — 1, or n, and it can prove very difficult to find generators of the vanishing ideal from 
them. For larger n it is also rather difficult to construct these polynomials. We constructed 
the following example. Let q be the simple Lie algebra of type A3. This Lie algebra has 
a 10- dimensional irreducible representation. Let xi, . . . , x 6 denote the matrices of the pos- 
itive root vectors of g in this representation, and let N be the space spanned by them. 
We tried to construct the above polynomials in Magma. But the system crashed when it 
exceeded the limit of 2GB of memory. The method based on elimination used 136 seconds 



4 



to construct defining polynomials (see Section [6]) . 



4 The semisimple case 

Let X G Ql{n, F) be a semisimple matrix. By G(X) we denote the smallest algebraic 
subgroup of GL(n, F) such that its Lie algebra contains X. In this section we consider the 
problem of obtaining defining polynomials for G(X). We will describe an algorithm based 
on the following theorem. 

Theorem 5 (Chevalley) Let a iy . . . ,a n G F be the eigenvalues of X . Let A G GL(n, F) 
be such that AX A' 1 = diag(ai, . . . , a n ). Denote this last matrix by Y . Set 

n 

A = {(ei,...,e n ) GZ" I ^0^ = 0}. 

Then 

G(Y) = {diag(ci, . . . ,c n ) \ Ci e F and c- 1 = 1 for all (e 1 ,...,e n ) G A}. 

i 

Furthermore, G(X) = A~ 1 G(Y)A. 

The first statement is [2], Chapter II, Proposition 2. The second statement is immediate. 
The first step in our algorithm will be to construct a finite extension F' D F containing 
the eigenvalues «j. Then by linear algebra we can construct a basis of the space 



A Q = {(ei, . . . , e n ) G Q n | ^ = °}- 



i=i 



From this we want to get a basis of A. First of all, by multiplying by suitable scalars we 
get a basis B of Aq whose elements have integral coordinates. By B we also denote the 
m x n-matrix having the elements of B as rows. Let S be the Smith normal form of B (cf. 
|12j). This means that we have unimodular matrices P, Q such that PBQ = S, with 

(d x •-• 0\ 

S = 

\ d m o ••• oy 

where the c?j are positive integers. 

Lemma 6 B is a basis of A z/ and on/y z/ a// are equal to 1 . 
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Proof. Suppose that all di = 1. Let A' C Z n be the subgroup generated by B. Then 
there is a surjective homomorphism / : Z n — > Z n_m with kernel A' ([12], Proposition 3.3 
of Chapter 8). Let v G A, and write d as a linear combination of elements of B with 
rational coefficients. After multiplying by a suitable integer s we see that sv G A'. But 
then = f(sv) = sf(v). Hence f(v) = and v G A'. We conclude that A = A'. 

Now suppose that B is a basis of A. We have a surjective homomorphism / : Z n — ► 
Z/diZ © ■ ■ • © Z/d m Z © Z n ~ m , with kernel A. Suppose that z is such that d { > 1. Then 
there is a f G Z n such that /(f) 7^ but f(div) =0. In other words, d(V G A. But then 
v G A, and we have a contradiction. □ 

Note that B = P~ 1 SQ~ 1 . This means that di times the z'-th row of Q _1 lies in the span 
of B, and hence in A. Therefore, the z-th row of Q' 1 itself lies in A, for 1 < % < m. Set 

(\ ••• o\ 

5" = 

v 1 ••• V 

and A = P~ 1 S'Q~ 1 . Then the rows of A belong to A, and form a basis of Aq. Moreover, 
the Smith normal form of A is 5". Hence by Lemma [61 the rows of A form a basis of A. 
We remark that there are efficient algorithms for computing the matrices P, Q (see |12j). 
Hence we can efficiently compute a basis of A, given a basis of Aq. 

Let Ai(X) be the associative algebra with one generated by X. Then Aj(X) is spanned 
by /, X, X 2 , . . . , X 1 , where t + 1 is the degree of the minimal polynomial of X. By |2J, §13, 
Theorem 10, G(X) C Aj(X). Let ol^Y be as in Theorem [51 We note that the minimal 
polynomial of a semisimple matrix is the square free part of its characteristic polynomial. 
In particular, it does not depend on the base field, but only on the coefficients of the 
matrix. So Aj(Y) is spanned by I, Y, ... , Y l . 

Let y = Y^i=Q$iY l G Aj(Y). Then by y(k,k) we denote the entry on position (k,k). 
By Theorem [5] we have that y G G(Y) if and only if y(k, k) ek = 1 for (ei, . . . , e n ) in 
a basis of A. Now set e' k = if > 0, and e' k = otherwise. Also e' k ' = e' k — e\~. Then 
jQ fc 2/(/c, k) £k = 1 if and only if Y\ k y{k, k) e k = Y\ k y{k, k) e k. In this we substitute y(k, k) = 

Sia l k . This yields a polynomial equation for the Si with coefficients in F' . By writing the 
coefficients as linear combinations of a basis of F' over F, we get polynomials p\, . . . ,p s , 
Pk G F[T , . . . , T t ], with the property that Yik vi^i ^Y k = 1 if and only if Pi(S , . . . , S t ) = 
for 1 < i < s. Let hi, . . . , h m be the totality of these polynomials that we get when we let 
(ei, . . . , e n ) run through a basis of A. Then y G G(Y) if and only if hi(So, ■ ■ ■ , S t ) = for 
1 < i < m. 

Note that y G G(Y) if and only if A~ 1 yA G G(X) (where A is as in Theorem [5]). But 
A~ x yA = J2i We conclude that 5^ G G(X) if and only if hi(S , . . . , S t ) = for 
1 < z < m. 

We summarise these findings in the following algorithm: 

Algorithm 7 Input: a semisimple matrix X G gl(n,F). 
Output: defining polynomials for G(X). 
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1. Construct a finite extension F' D F containing the eigenvalues of X . 

2. Compute a basis of A. 

3. Construct polynomials hi, . . . , h m G F[T , . . . ,T t ] with the property that £V G 
G(X) if and only if hi(5o, . . . , 5 t ) — for 1 < i < m. 

4. Set M = Y?i=v T i xi and let a ij e F be such that T k = £V a^M(i,j). 

5. Let Xij for 1 < i, j < n be indeterminates , and consider the substitution T k \— > 
Y^ij a ij x ij- Let ip : F[T , . . . , T t ] — > fre the corresponding ring homomorphism. 

6. Let g 1 ,...,g r G i* 1 ^-] be linear polynomials with the property that a G Ql(n,F) lies 
in Aj{X) if and only if gi(a) = for 1 < i < r. 

7. Return {g ± , . . . ,g r } U {ip{hi), . . .,ip(h m )}. 

Proof. Steps 1, 2, and 3 have already been commented on. Since the X 1 for < i < t 
are linearly independent, the coefficients afj in Step 4 exist (but they are not necessarily 
unique). Let a = (aij) G GL(n, F) be such that gi(a) = for all i. Then a G Aj(X) so 
a = Ylii^iX 1 - By substituting T k \— > 5^ in the equation T k = . a^-M(i, j) we see that 
<5 fc = £V . a^Oij. It follows that a G G(X) if and only if 

h AYl a % ai ^ ■ ■ ■ ' XI a <i a «) = 

for 1 < r < m. But this is equivalent to (p(h r )(a) = for 1 < r < m. We conclude that 
this algorithm returns defining polynomials for G(X). □ 

Example 8 Let 



Then Aj(X) is spanned by I,X. The eigenvalues of X are %,—%, and hence A is spanned 
by (1, 1). Let Y = diag(i, —i) and y = 5 I + 5iY, then y(l, 1) = 5 + <M, 2/(2, 2) = 5 — 
Since dim A = 1, there is only one equation that needs to be satisfied, i.e., y(l, l)y(2, 2) = 1, 
or ^o + ^i = !• Now 

M = T I + T X X = 

So we get the substitution T = j u and T x = rr 12 . The equation above then transforms 
to x\y + x\ 2 — 1 = 0. Furthermore, the linear equations of Step 6 of the algorithm are 

^21 + %V2 = 0, X n - X 2 2 = 0. 
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5 An application 



Let G C GL(n, F) be an algebraic group, defined over F. Then there is a unique maximal 
unipotent subgroup, u(G), called the unipotent radical of G. Furthermore, there is a 
reductive subgroup H C G with G = H x u(G). In [7] algorithms are given for finding 
defining polynomials for if and u(G). However, these algorithms do not seem to be very 
practical. In this section we describe a method based on the algorithms of the previous 
sections. 

Set q = Lie(Cr). Let s denote the solvable radical of q, and I a Levi subalgebra of g, 
and n the largest ideal of q consisting of nilpotent elements. Then by [T|, Chapter V, §4, 
Proposition 5 (see also [H]), s has a commutative subalgebra consisting of semisimple 
elements, with the following properties 

1. s = d + n (semidirect sum), 

2. M] = o, 

3. n is the set consisting of all nilpotent elements of s. 

Now H and u(G) are the algebraic groups corresponding to [ + U and n respectively. 
There are algorithms to compute 1 and s (cf. |5j). In the remainder of this section we 
describe how to compute and n. Then by using the algorithms in the first part of this 
paper, we can compute defining polynomials for H and u(G). 

Lemma 9 Set I) = C S (V), the centralizer of d in s. Then I) is a Cartan subalgebra of 5, 
andt) = Q + C n (X>). 

Proof. Let x G C s (p), then we can write x = y + z, where y G and z G n. So for 
any u G 9 we get = [u,x] = [u,y] + [u,z] = [u,z]. It follows that z G C s (d). Hence 
f) = + C n (d). Since n is nilpotent and [0, C n (())] = it also follows that f) is nilpotent. 

Now let x G 5. Then x is a semisimple linear transformation. This implies that also 
&dx : s — » 5 is semisimple. Hence C s (0) = So(adc)), where 

So(adO) = {y G s | (ada;) dims (?/) = for all x G d}. 

By [13], Theorem 4.4.4.8 every Cartan subalgebra of So(adO) is a Cartan subalgebra of s. 
But as seen above f) = So(adO) is nilpotent. So it is its own Cartan subalgebra. □ 

Lemma 10 Let h' be a Cartan subalgebra of s. Then there is an abelian subalgebra d' C s 
consisting of semisimple elements, such that 5 = O' + n (semidirect sum), and I)' = C B (d') = 

Proof. Let D be the group of automorphisms of s generated by exp(adx) for x G [s,s]. Set 
f) = C s (d), which is a Cartan subalgebra of s by Lemma [91 By [I], Chapter VI, Proposition 
19, there is a a G D such that cr(fj) = h'. Set c)' = o~(d). In order to prove that d' consists 
of semisimple elements we may assume that o = exp(ada;). Then a(u) = (exp x)u(exp x)" 1 
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(by [9J. Chapter IX. (38). note that [s.s] C n. so cxp(.r) is a well defined endomorphism). So 
d' is obtained from d by conjugation by a fixed element, and hence it consists of semisimple 
elements. The other properties follow from the analogous properties of f) and the fact that 
a is an automorphism. □ 
Let r be a finite-dimensional Lie algebra, and let a C t be a nilpotent subalgebra. Then 
r has a Fitting decomposition with respect to the adjoint action of a (cf. [9], [5]). This 
decomposition is written r = r (a) © ti(a). These are called the Fitting 0-component and 
Fitting 1-component respectively. We have ti(a) = flj >0 [a*, t], where by [a\ r] we denote the 
space [a, [a, • • • , [a, r] • • •]] (i factors a). Based on this there is a straightforward algorithm 
for computing a basis of ti(a) (see [5]). 

Lemma 11 Let h' and d' be as in the previous lemma. Let a±, . . . , a r be a basis of h' and 
let ai = Si + Hi be the Jordan decomposition of aj. Then d' is spanned by s%, . . . ,s r . Let 
£i(f)') be the Fitting 1-component of s with respect to the adjoint action of h'. Then n is 
spanned by the along with Si(h'). 

Proof. Let h G f)'; then by Lemma [10], there are s£i)' and n G C n (p') with h = s + n. 
But then s is semisimple, n is nilpotent and [s, n] = 0. It follows that h = s + n is the 
Jordan decomposition of h. In particular, all lie in d' and all lie in C n (0'). Also 

h = J2i a i a i = (Y,i a i s i) + {J2i a i n i)- Tnis implies that J2i a i s i = s ancl J2i a i n i = n - 
Therefore, the Sj span 0', and the span C n (D'). Note that So(f)') = f}' as h' is a Cartan 
subalgebra of s (cf. [9], Chapter III, Proposition 1). Also Si(h') C [s,s] C n. From this we 
get the last statement. □ 
We note that there are efficient algorithms for computing a Cartan subalgebra of a 
Lie algebra of characteristic (cf. |5J). So the previous lemma yields a straightforward 
algorithm for computing bases of D' and n. 



6 Practical experiences 

In this section we report on practical experiences with an implementation of the algo- 
rithms in the computer algebra system Magma. The computations were done on a 2GHz 
processor, with 2GB of memory. 

First we consider the algorithm from Section [3j This algorithm takes as input a set 
of nilpotent matrices that span a Lie algebra, and returns defining polynomials for the 
corresponding algebraic group. In order to generate input we have constructed irreducible 
representations of several simple Lie algebras, and taken the matrices of the positive root 
vectors. The running times of the algorithm are listed in Table [U From the table we see 
that the algorithm performs quite well for small matrices. However, the running times 
increase sharply when the size of the matrices increases. This is due to the Grobner basis 
calculation, used for the elimination of variables. 

Next we consider the algorithm of Section HI Here the input is a single semisimple 
matrix. As input we have used the companion matrices of several square free polynomials 
(except the first one these are taken from the database in [ID])- The running times are 
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type 


dim rep. 


dim iV 


time G(N) 


B 2 


5 


4 


0.08 


B 2 


10 


4 


43 


G 2 


7 


6 


3.45 


A 3 


10 


6 


136 


c 3 


14 


9 


oo 


D 4 


8 


12 


oo 



Table 1: Running times for the algorithm of Section [3] The first column lists the type of 
the simple Lie algebra, the second column has the dimension of the irreducible module, 
the third column the dimension of the Lie algebra spanned by the positive root vectors. 
The last column displays the time (in seconds) for constructing defining polynomials for 
the group. An oo means that the algorithm did not complete with 2GB of memory. 



polynomial 


dimension splitting field 


time G(X) 


x 8 - x 7 + x 6 + 2x 5 - 3x 4 + 4x 3 + 2 


360 


83 


x 6 - 2x A + x 2 - 2x - 1 


360 


293 


x 8 + 12x 6 + 50x 4 + 83x 2 + 43 


384 


89 


x 8 - 2x 6 + 7x 4 - 8a; 2 - Ax + 7 


576 


1463 



Table 2: Running times for the algorithm of Section HI The first column lists the polyno- 
mial, and the second column the degree of its splitting field. The last column displays the 
time (in seconds) for constructing defining polynomials for the group corresponding to the 
companion matrix of the polynomial. 

displayed in Table [2j As was to be expected the algorithm spends a lot more time if the 
degree of the splitting field gets larger. 

Finally we consider the algorithm of Section [2J For this algorithm we have constructed 
several rather different inputs, which we describe separately. 

Let A be an algebra, and Der(A) the Lie algebra of its derivations. By [2], Chapter II, 
§14, Theorem 16, we know that Der(A) is the Lie algebra of the automorphism group of 
A. Here we consider the algebra A of quaternions over Q. This algebra has basis k, 
where 1 is the identity. The other product relations are ij = k, ji = —k, ik = —j, ki = j, 
jk = i, kj = —i. Then Der(A) is spanned by 





/o 








o\ 




A) 








o\ 




/o 





o\ 










1 
















1 













Xi = 





-1 






























1 




V' 












v° 


-1 











-1 


0/ 



(where we use the column convention to represent a linear map with respect to the basis 
1, i, j, k of A). Let Hi be the smallest algebraic group containing Xj. The x^ are semisimple, 



10 



so equations for the Hi can be computed using the algorithm of Section HI Let G be the 
smallest algebraic subgroup of GL(4, C) containing Hi, H 2 . The Lie algebra of G contains 
Xi,x%- But then it contains also x%. It follows that its Lie algebra is equal to Der(A). 
Hence G is the connected component of the identity of the automorphism group of A. The 
computation of equations for G (given those of Hi and H 2 ) cost 8.3 seconds. 
The next example has the Lie algebra consisting of the matrices 



/* 


* * 




* 


* * 


* 





* 


* 




* 





Let Hi denote the group corresponding to the subalgebra consisting of strictly upper 
triangular matrices. By Hij we denote the group corresponding to the subalgebra spanned 
by ey (where denotes the matrix with 1 on position and zeros elsewhere). Let 
Gi denote the smallest algebraic group containing both Hi and H 2> i- The construction 
of defining polynomials for Gi cost 0.99 seconds. Let G 2 denote the smallest algebraic 
group containing Gi and H^. The computation of defining polynomials for G 2 cost 2.4 
seconds. Let G 3 be the smallest algebraic group containing G 2 and H11. The computation 
of defining polynomials for G3 cost 5.4 seconds. Let G = G4 denote the smallest algebraic 
group containing G3 and -£^3,3. The computation of defining polynomials for G\ cost 15 
seconds. We note that G is equal to the algebraic group corresponding to the whole Lie 
algebra. 

Let L be the Lie algebra of type B 2 in its 4-dimensional representation. Let Hi and H 2 
denote the groups corresponding to the subalgebra spanned by the positive and negative 
root vectors respectively. Let G be the smallest algebraic group containing both Hi and 
H 2 . The computation of defining polynomials for G cost 173 seconds. 

Finally, let L be the simple Lie algebra of type A% in its 4-dimensional representation. 
Let Hi, H 2 and G be as in the previous example. The algorithm did not manage to compute 
defining polynomials for G within 2GB. 

Concluding we can say that the algorithms do work for small examples. However, 
due to the Grobner basis computations, needing a lot of memory and running time, the 
algorithms are not yet practical for larger examples. 
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